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Abstract 



We re-analyze the Schwinger boson mean field theory (SBMFT) for 
Heisenberg spin models on the cubic lattice. We find that the second order 
phase transition point for magnetic ordering previously reported corresponds 
to a local maximum of the free energy functional. For both ferromagnetic and 
antiferromagnetic Heisenberg models with spin S > Sc, where Sc < 1/2, the 
mean field transitions are first order from the magnetically long-ranged or- 
dered phase to the completely uncorrelated phase. In addition to erroneously 
giving a first order transition for magnetic ordering, the mean field theory 
does not include a phase with finite short-range correlation, thus negating 
one of the prime advantages of SBMFT. The relevance of these pathologies 
to other situations beyond the cubic lattice is discussed. 
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I. INTRODUCTION 



Schwinger boson representation is by now a well known tool in the study of quantum spin 
systems In this method, each spin operator is represented by two Schwinger bosons. 

In contrast to Holstein-Primakoff bosons [Q, the Schwinger boson representation of spin 
operators does not involve square roots. Furthermore, while the vacuum state of Holstein- 
Primakoff bosons has to be a broken symmetry state, this is not the case for Schwinger 
bosons. Instead, long-ranged order (LRO) is a result of condensation of the Schwinger 
bosons [^0. Thus, using Schwinger bosons allows one to naturally address both the mag- 
netically disordered and ordered states. 

Using the Schwinger boson representation, the quantum spin Hamiltonian is mapped 
onto an interacting boson Hamiltonian, which can then be approximated by a self-consistent 
non- interacting Hamiltonian using mean field theory (MFT). In such a mean field theory, a 
non-zero mean field amplitude corresponds to short-ranged order (SRO). Long range order 
exists only if for a given temperature, the mean field Hamiltonian gives rise to Bose conden- 
sation. The ability of the mean field theory to address both short-ranged and long-ranged 
order is perhaps its most appealing feature, especially for ID and 2D systems, where exact 
solutions has shown that LRO is destroyed at any finite temperature. For the Heisenberg 
ferromagnetic or antiferromagnetic Hamiltonian on a d-dimensional hypercubic lattice, the 
results reported in literature are briefly summarized as follows. In both ID and 2D, in agree- 
ment with exact results, there is no LRO at any finite temperature. However, SRO persists 
until a certain temperature is exceeded In 3D, LRO is stable at low temperature, 

and there is a continuous phase transition into the SRO state at some critical temperature. 
Within Schwinger boson MFT (SBMFT), the order- disorder transition has the usual mean 
field exponents [0. 

In this paper, we re-examine the SBMFT for the Heisenberg spin Hamiltonian on the 3D 
cubic lattice. We find that, provided the spin S is greater than some critical value Sc {< 1/2), 
the previously reported mean field solution that exhibits a continuous phase transition for 
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magnetic LRO is actually an unstable solution corresponding to a local maximum of the free 
energy. Instead, the minimum free energy solution shows a first order transition. Since it 
is well established that the Heisenberg Hamiltonian has a continuous transition in 3D, our 
analysis raises serious question about the applicability of the SBMFT in 3D. Furthermore, 
the first order transition is into a state with zero mean field amplitude, i.e., without any 
correlation. Thus, the supposed advantage of SBMFT of being able to describe SRO does 
not actually exist. These pathologies of the SBMFT are not restricted to 3D, but should 
appear whenever the LRO can be destroyed by tuning some parameter. An example of which 
is the bilayered Heisenberg antiferromagnet at T = where as the inter-layer coupling is 
increased, SBMFT also shows a first order transition into the uncorrelated state 0. 

The rest of this paper is organized as follows. In Section II and III, we will discuss the 
SBMFT solutions for the Heisenberg ferromagnet and antiferromagnet separately, showing 
that the same pathologies are seen in both cases. These results will be summarized in Section 
IV, where we will also discuss their relevance to the bilayered Heisenberg antiferromagnet. 



II. FERROMAGNETIC HEISENBERG MODEL 

In this section we discuss the SBMFT for the FM Heisenberg model. We first briefly 
review the Schwinger boson representation and its mean field equations |T]-§[. We then 
discuss the mean field solutions of the model with emphasis on 3D. 

We study a spin-S Heisenberg ferromagnet in a d-dimensional hypercubic lattice with 
periodic boundary conditions, 

H=-jJ2S^-Sv (1) 

<i,j> 

where the exchange constant J > and the sum is over all the nearest-neighbor pairs on the 
lattice. In the Schwinger boson representation, a spin operator is represented by two types 
of bosons, 

-^r = ^kV^t (2) 
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S^^a^^aii, =aj^ait, (3) 

These bosons should satisfy the constraints in their Fock space as follows, ctla^if^ ~ 25", 
with 25" = 1, 2, ... for physical systems. By introducing a bond operator, Bjj — | a\^ajcr, 
the SU(2) FM Heisenberg Hamiltonian can be written as, 

Hfm-b --'^JYI - 4^'^- • +^-5' + E - 2-5)- (4) 

<ij> i c 

In the above Eqn., :: denotes a normal ordering of the bosonic operators, N is the 
number of sites in the lattice, and z = 2 d is the nearest neighbor coordinate. A Lagrangian 
multipher field Aj for the constraint on site % is introduced. Note that the FM Heisenberg 
model is now transformed to an interacting boson model. 

To proceed further, a bond mean field amplitude B =< Bf- >=< Bij >, assumed to be 
real and uniform, is introduced. Neglecting fluctuations, the following non-interacting mean 
field Hamiltonian is obtained., 

Hmf = A alui, - 2S) - 2 J J] B{Bj^ + S,,) + + NzJB\ (5) 

i c ij 

In Hmf the local constraint at every site is replaced by a global one, and a non-zero value of B 
indicates short-range spin-spin correlation. The mean field Hamiltonian can be diagonalized, 
and we obtain 

Hmf = + NzJB' - 2SXN + ^^^L^a?.' (6) 

where the energy dispersion is given by a;^ = JzB{ej^ + A), with eg = ^ — e^^-^), and 

The free energy per site is then given by 

^ ^j^2 _2s\-A.y iog( — L_), (7) 

k 

where P — -j^, with ks the Boltzmann constant and T the temperature. The amplitudes 
of the mean fields B and A are obtained by minimizing the free energy, which leads to two 
coupled mean field eqns.. 



where (3 = l/fc^T, and T = T/(JzB) acts as an effective temperature. In these eqns., —A 
may be viewed as an effective chemical potential. Since > 0, we have A > 0. Magnetic 
long range order corresponds to Schwinger boson condensation at A; = Since bosons 

do not condense at finite temperature in ID or 2D, the ID model and the 2D model are 
disordered at any temperature. In three and higher dimensions, one expects a finite critical 
temperature Tc, below which SBMFT gives a magnetic ordered state. 

To study the possible magnetically ordered phase, the Schwinger boson condensate den- 
sity p =< a^_Q^ctfc=oo- > /-^ is introduced, and the coupled mean field eqns. are rewritten 
as. 

In the Bose condensed state, p > 0, A = 0, and i? 7^ 0; In the disordered state, p = 0, and 
A > 0. Sarker et al. determined the critical temperature by requiring A = and p = 0. 
This transition point will be labeled as "a", with the effective temperature T^. By expanding 
equations 10 and 11 about T^, the transition was determined to be continuous with mean 
field exponents. Upon re-examination of the mean field equations, we have, however, found 
the transition point to be physically inaccessible if the spin S is greater than some Sc- To 
understand this, we note that these equations are well-behaved in term of T. For example, 
the quantity B is smooth and monotomically decreasing as T increases as one would expect. 
This is shown in fig. 1 for the case of S* = 1/2. At T = 0, -B is finite, and A = 0. As T 
increases, B decreases monotonously. For T > T^, the state is disordered with p = 0; And 



at T < Ta, the state is ordered with p > 0. On the other hand, the physical temperature 
T = JzBT, and since B decreases with increasing with T, it is possible for T to decrease as 
T increases. This is shown in fig. 2. This would mean for example that B would increase 
with increasing T, which is physically nonsensical. If the second order transition point 
belongs to such an unphysical T, it will not be accessible. 

In ref.2, it was implicitly assumed that the mean field eqns. 9 and 10 have a two branch 
structure, so that at high T, S = is the only solution, but below some T', a second solution 
with B ^ appears. Instead, mean field eqns. actually have three branches of solutions for 
5 or T as functions of T. One is the trivial solution (upper branch), 5 = 0, or T ^ oo at all 
T. The second branch solutions (lower branch) exist for all T < Tf, with T monotomically 
increasing with T. For Tg < T < Tt, there is yet a third branch of solutions (middle branch), 
where T goes to oo at T = Tg and decreases as T increases until it merges and terminates 
with the second branch at T = T^. 

The solutions to the mean field equations correspond to extrema of the free energy. For 
T < Tf, there are more than one solutions, and the one giving the lowest free energy should 
be chosen. From the structure of the T vs. T curve, we can deduce that as T is reduced, the 
free energy functional should behave as follows. For T > Tj, there is only one extremum, 
which is at i? = 0. Between Tt and Tg, where there are three solutions, there is a local 
minimum at B — 0, a local maximum at B — B^ (middle branch) and a secondary local 
minimum at B — B2 (lower branch), with B2 > B3. As T decreases in the regime, S3 moves 
towards and at T = T^, the curvature at B = changes sign and becomes a local maximum 
below Tg^ leaving B2 as the only minimum. Given this, at some temperature Tf, between Tg 
and Tt, the free energy at B2 will become lower than at B — 0. Thus, at T^^ B will take a 
jump from to a finite value of B2. In Figs. 3-5, we show plots of free energy / vs. B to 
illustrate the above. In these plots, the f{B) is calculated by solving eq. 10 for A with B 
as a variable and then evaluating / in eq. 7. 

Thus, we see that in general, SBMFT gives a first order transition in B. Since i? is a 
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measure of the short range order, the transition in B is fictitious and whether it is first 
order or continuous is not of any real significance. More importantly is where the Bose 
condensation transition occurs relative to this transition. Since in terms of T, eq. 10 is just 
that for Bose- Einstein condensation, ~ 5'^/'' for d > 2, as 5" ^ 0, and hence so does 
the physical temperature T^. Thus for small S, Bose condensation occurs at a temperature 
where the lowest branch is the stable solution, and the continuous transition at accessible 
within SBMFT. As S increases, however, Ta continues to increase and for S > Sc ^ 0.15, 
Ta > Tb, where the solution at the lower branch becomes meta-stable. At even larger S (e.g. 
S = 1/2), Ta becomes > Tj, and the solution becomes unstable. Either way, for S > Sc, as 
the temperature increases from zero, the system undergoes a discontinuous transition at 
from the LRO state with finite condensate density to the uncorrelated state with B — 0. 
There are two problems with this. The first is that result of a first order transition in 
magnetic long range order for the Heisenberg model in 3D is wrong. The second is that one 
of the major advantages of SBMFT is supposedly its applicability to both LRO and SRO 
state. But in fact, the stable solution is either completely uncorrelated or has LRO, just like 
standard Bethe- Weiss MET. 

III. ANTIFERROMAGNETIC HEISENBERG MODEL 

In this section, we discuss the SBMET for the AEM Heisenberg model described by the 
Hamiltonian, 



where J > 0, and all other notations are the same as in the previous section for the EM 
case. We find that the same pathologies in the SBMET solution for EM model are present 
also in the AEM model. 

The hypercubic lattice is divided into two sub- lattices A and 5. In a nearest neighbor 
pair < ij >, ieA and jEB. As in the EM model, two Schwinger bosons Ojo- with a — ±1, 




(12) 



<i,3> 



7 



at each site i, and a bond field = | X^o- '^^L^j-o- introduced. Aij is anti-symmetric 
with respect to the interchange of i and j. A spin rotation by tt about the y-axis is made 
at all the sites j on the sub-lattice B, so that aj^ — >aji, aji — ^Oj-f. The bond operator is 
then transformed as Ajj — > Ajj — | at^aj-^. This canonical transformation preserves 

the constraint aL'^«cr = 25" = 0, 1, 2, 

In terms of bosonic operators, the AFM Heisenberg model reads, 

— E(4^^^-2'5')- (13) 

<ij> 

The mean field equations of the AFM model may be derived similarly to the FM model. 
A global Lagrangian multiplier A to replace the constraint on every site, and a real and 
uniform mean field amplitude < A^j >—< Ajj >— A are introduced. A describes short- 
range AFM correlations. After some algebra, the mean field Hamiltonian can be written 
as: 

Hmf = Eo + XJ2 aL«fe^ Y~ 5Z ^''^"'l^^-ka + ^-c), (14) 

with 7fe = i e'^'^^ = 1 - ^fe' and Eo = ^NzJS^ - 2XNS + JA^Nz. 

Hmf is diagonalized by using the standard Bogoliubov transformation, yielding 

Hmf ^Eo-XN + Y,M4<^k + PlPk + 1), (15) 

k 

where lOj^ — (A^ — {JzA'j^yy^^. To find the optimal values of the mean fields, the free 
energy is differentiated with respect to A and A, and the SBMFT equations for the AFM 
model are obtained, 

+ 2 " TV ^ (At2 - 72)1/2 (g^(/.2-^2)l/2 _ + 2) ' (^^^ 

AT ^ (/X2 _ ^2) 1/2 V g/3(;.2_^2)i/2 _ ^ ^ 2 y ' ^ ^ 

where we have defined an effective inverse temperature /3 = JAzj3 and an effective 
chemical potential ji — -j^. To describe the possible Schwinger boson condensed state. 



condensation density p is introduced similar to that in the FM model, and the summations 
are replaced by integrals. The mean field eqns. become, 



„ 1 f d'^k 1 

S + 7; = P 



(e/3(M'-7p^/' -1 ^ 2)' ^^^^ 



2 ^ y (27r)«'(/x2-7|)V2Ve/3(M^-7p 



We now discuss the solutions of these eqns. Note that 7H<1, so we have ii>l. Similar 
to the FM case, p = and // > 1 in the disordered phase, and p > and // = 1 in the 
Schwinger boson condensation state, which corresponds to the Neel ordering phase of the 
spin system. The results are very similar to those for the 3D FM model. 

The solutions for the mean field equations are smooth and monotonic functions of the 
effective temperature T, but not of the physical temperature T. In Fig. 6, we show the 
calculated T as a function of T from the mean field eqns. for a typical spin S = 1/2. There 
are three branches of solutions in T. The A — Q branch, and a middle and a lower branches. 
The Bose condensation point "a" (where ix — 1 and p = 0) is again on the upper branch 
corresponding to a maximum in the free energy. We have calculated the free energy and 
found that the stable mean field solution jumps from a high T phase in the A — Q branch 
to a low T phase in the lower branch with a finite condensation density at a temperature 
T — Th < Ta- Therefore, the transition for LRO is first order. The qualitative feature of the 
free energy as a function of A is similar to those in Fig. 3-5 for the FM case. The critical 
spin for which Ta — Tj, for the AFM model is found to be Sq — 0.32, larger than that in 
the FM model, but still substantially smaller than 1/2. For S > Sc the ordering transition 
within SBMFT is first order and for S < Sc, the transition is second order. 



IV. CONCLUSIONS 



We have re-examined the Schwinger boson MFT for quantum spin systems and found 
that for both the FM and AFM Heisenberg Hamiltonian, the mean field equations possess 
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in general three branches of solutions in the intermediate temperature regime. As a conse- 
quence, the system in 3D undergoes a first order transition from the ordered state to the 
completely uncorrelated state unless the spin S is very small. In addition to erroneously 
predicting a first order transition in magnetic ordering, there is no stable solution with fi- 
nite SRO without LRO, the existence of which is supposed to be one of the advantages of 
SBMFT over Bethe- Weiss MFT or spin wave theory. 

The three branch structure of the mean field solutions is not restricted to 3D. It is 
the origin of the reported discontinuous jump in the short-range MF amplitude A for the 
Heisenberg AFM on 2D square lattice as the temperature is increased . For the 2D square 
lattice, however, the effective temperature for Rose condensation = 0, which means 
Ta < Tb for all S. As a result, SBMFT is capable in that case of describing the SRO state 
at low temperature. It should not be assumed though that the pathologies associated with 
SBMFT reported in this paper is only relevant for d > 2. In fact, it will apply whenever 
the LRO state is destroyed by tuning some parameter if that parameter is sufficiently large. 
One example of this is the bilayered Heisenberg AFM Hamiltonian in 2D, which received 
interest a few years ago due to its relevance to YBCO In the bilayered system, the 

in-plane AF magnetic LRO is destroyed by sufficiently large inter-layer coupling which favors 
having inter-plane neighbors forming singlets P,|6|,PHl^. Using SBMFT, Ng et al ^ found 
for S > Sc — 0.35, the transition is first order from finite condensate density to the state of 
disconnected inter-plane singlets with zero in-plane correlations. This result contradicted the 
hypothesis that the 2D bilayered AFM at T = is in the same universality class as the 3D 
non-linear sigma model at finite T, a hypothesis partially confirmed by finite sized numerical 



calculations for S* = 1/2 ||14|. One possibility is that SBMFT is qualitatively correct, but 
that the value of Sc calculated is too small. Subsequently, Gelfand et al provided further 
evidence that the transition for arbitrary 5* is continuous and that the SBMFT prediction 
of first order transition is incorrect. However, it was not clear why SBMFT gave such a 
qualitatively incorrect result. From the present work, we can now understand that the first 
order transition reported there is an artifact of the same pathologies of the SBMFT reported 
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in this paper. 

Finally, while we have focused on SU(2) Heisenberg Hamiltonian in this paper, the results 
are a consequence of the structure of the mean field equations, and therefore will be valid 
for SU(N) generalization of the Heisenberg Hamiltonian for which the SBMFT approach 
above can be applied. 
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FIG. 1. Effective temperature (T) dependance of B for three dimensional spin 1/2 Heisenberg 
ferromagnetic Hamiltonian. Bose condensation occurs at Ta- 
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FIG. 2. Effective temperature (T) as a function of physical temperature {T/J) for three di- 
mensional spin 1/2 Heisenberg ferromagnetic Hamiltonian. 



13 



0.5 1.0 



B 

FIG. 3. SB mean field free energy at T = 1.3333 J > Tf for three dimensional spin 1/2 
Heisenberg ferromagnetic Hamiltonian. 
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FIG. 4. SB mean field free energy at < T = 0.9753 J < Tt for three dimensional spin 1/2 
Heisenberg ferromagnetic Hamiltonian. 
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FIG. 5. SB mean field free energy at T = 0.8333 J < Tj, for three dimensional spin 1/2 Heisen- 
berg ferromagnetic Hamiltonian. 
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FIG. 6. Effective temperature (T) as a function of physical temperature (T/J) for three di- 
mensional spin 1/2 Heisenberg antiferromagnetic Hamiltonian. 
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